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Abstract 

We extract the high-twist contribution to the neutrino-nucleon structure function 
xF^ u+u ^ N from the analysis of the data collected by the IHEP- JINR Neutrino Detector 
in the runs with the focused neutrino beams at the IHEP 70 GeV proton synchrotron. 
The analysis is performed within the infrared renormalon (IRR) model of high twists 
in order to extract the normalization parameter of the model. From the NLO QCD 
fit to our data we obtained the value of the IRR model normalization parameter A| = 
0.69 ± 0.37 (exp) ± 0.16 (theor) GeV 2 . We also obtained A§ = 0.36 ± 0.22 (exp) ± 
0.12 (theor) GeV 2 from a similar fit to the CCFR data. The average of both results is 
A§ = 0.44 ±0.19 (exp) GeV 2 . 

PACS numbers: 13.60.Hb,06.20. Jr,12.38.Bx 
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1. Attempts to extract the high twist (HT) contributions to the neutrino-nucleon deep- 
inelastic scattering (DIS) structure functions started many years ago |J, but have not lead 
to the ultimate answer up to now. The main difficulty in this study is that due to the 
linear rise of the total interaction cross-section with the incident neutrino energy E u , the 
largest data samples have been collected in experiments at relatively high neutrino energies 
E v > 50 GeV. The region of the small momentum transfered Q, which is most relevant 
for the study of the HT effects, is rather poorly populated by the data points coming from 
these experiments because of kinematical and/or methodical cuts. At lower neutrino energies 
experiments with very high luminosity are necessary to achieve the statistical precision in 
the structure function measurements sufficient for the quantitative estimation of the HT 
contribution. 

The first results on the twist-4 contribution to the neutrino-nucleon structure function 
xF^ N extracted from the analysis of data collected in a single experiment were reported 
in Ref. |2|]. The indication on a negative sign of this contribution given in this paper was 
later confirmed with a better precision in Ref. ||. Nevertheless, the experimental errors 
were large in both cases which did not allow for a conclusive comparison with the available 
theoretical models of HT. Later the CCFR experiment at Fermilab collected large statistical 
data sample Q, which allowed for more precise determination of the twist-4 contribution to 
xF% N 0, m, but the precision is still poor. 

The data of Ref. [Q from the IHEP-JINR Neutrino Detector can be used to improve our 
knowledge of the HT contribution to the DIS structure functions. This experiment used a 
neutrino beam of relatively low energy [E v < 30 GeV), but collected rather large statis- 
tics (5987 neutrino and 741 antineutrino charged- current (CC) interactions). The lowest 
Q 2 is 0.55 GeV 2 and the HT contribution would clearly manifest itself as a power-like cor- 
rection to the logarithmic-like leading twist (LT) dependence of the structure functions on 
Q. Meanwhile the Q range spanned by the data is limited (maximal Q 2 is 20 GeV 2 ) and 
for this reason the simultaneous determination of the power-like and logarithmic-like terms 
is difficult. In the analysis of Ref. |7j we fixed the HT contribution as it was defined from 
other experiments and performed the next-to-leading-order (NLO) QCD analysis of our data 
in order to constraint the LT contribution. The value of the strong coupling constant a s , 
which mainly governs the Q dependence of the LT term, was determined from this analysis 
as a(M z ) = 0.123lg;{$j. 

At the same time the Q dependence of the LT contribution can be well constrained using 
the data from other experiments. The world average of a s (Mz) is known with the precision 
of about 0.003 || and one can perform the analysis complimentary to the one of Ref. |7]], i.e., 
fix the value of a s (Mz) at the world average and try to extract the HT contribution from the 
data. An additional error estimated as variation of the results of the fit under variation of 
the world average of a s (Mz) within its uncertainty should be ascribed thereafter. Meanwhile 
in most cases the total error in the HT contribution extracted using this approach would be 
less as compared with the results of the simultaneous fit of the HT and the LT terms to the 
data. The reduction of the error depends on the ratio of the error in a s (Mz) obtained in 
such simultaneous fit to the error in the world average. For our data this ratio is larger than 
3 and for this reason we extracted the HT contribution value from these data using the fit 
with a fixed value of a s . 

2. The analysis is based on the data collected with three independent exposures of the 
IHEP-JINR Neutrino Detector |J to the wide-band neutrino and antineutrino beams [TIJ 
of Serpukhov U-70 accelerator. The exposure to the antineutrino beam (77^-exposure) was 
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performed at the proton beam energy E p = 70 GeV, whereas the two z/^-exposures were 
carried out at E p = 70 GeV and at E p = 67 GeV. The energy of the selected (V^) CC 
events was in the range of 6 < E V M) < 28 GeV. The experimental set-up and the selection 



criteria of CC neutrino and antineutrino interactions are discussed in Ref. [(LTJ]. The F 2 + " ]S 



and xF^ structure functions of nucleon have been measured as a function of x averaged 
over all Q 2 permissible (the details of experimental procedures are described in Ref. 0). 

These data were analyzed in the NLO QCD approximation in the modified minimal- 
subtraction (MS) renormalization-factorization scheme. The partons evolution code applied 
in this analysis was used earlier for the global fit of the parton distribution functions ||12|| . 
The boundary parton distributions were chosen in the form 

x Pm {x, Q ) = A m x aNS {l - x) feNS , x Ps {x, Q ) = A s (l - x) bs , 

xp G (x,Q )=A G (l-x) bG (1) 

at Ql = 0.5 GeV 2 , where indices NS, S, and G correspond to non-singlet, singlet, and gluon 
distributions, respectively. These distributions were substituted in the expressions for the 
LT contributions to F 2 $ 

xFt +m ' LT (x, Q) = f -Ci(z, Q)p NS (x/z, Q), 

J X z 

p (u W)N,LT^ Q) = f 1 dz {cq ^ Q) [pNs{x/ ^ Q) + ps{x/ ^ Q)] + 
J x z 

+C G (z,Q)p G {x/z,Q)} (2) 

where C(z, Q) are the perturbative QCD coefficient functions in the MS scheme. The pa- 
rameter Ans was calculated using the constraint dxq^s = 3, and the parameter A G - from 
the momentum-conservation constraint, while other parameters of Eq. (|l]) were fitted to the 
data. The form of Eq. ([]]) was checked to be flexible enough, i.e., its complication did not 
lead to the improvement of the fit. 

The target mass (TM) corrections of 0(M 2 /Q 2 ) } as they are given in Ref. [01] 



were 



applied to the LT contribution. The HT contribution was parameterized in the additive form 

and within the infrared renormalon model (IRR) [14]]. In this model the HT contribution is 

connected with the LT one by the known coefficient function and the only free parameter 

of the model is related to the total normalization. In particular, the HT contribution to 
xF i»+W reads |]| jig 



H 3 (x, Q) = A' 2 {F» N ) / -Cl^Wpssix/z, Q), (3) 

J x z 

where C^ RR is the IRR model coefficient function and A' 2 (F^ N ) defines the total normal- 
ization. The HT contribution to E^ contains the non-singlet term similar to Eq.([3]) 
with the normalization parameter A' 2 (F2 N ) and the respective coefficient function C 2 RR . In 



addition, the singlet and gluon terms calculated in Ref. |L7j also come to the expression for 



the HT contribution to F 2 as it is given by the IRR model, but these terms are relevant for 
small x only and for this reason we used the non-singlet approximation for the calculation 

7 b ' 



of the IRR contribution to E 2 as well as for xE^ u+u ^ N . Following Ref. W71, we describe 
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the general normalization of the HT contributions to 7^,3 by the parameters A2,3, which are 
connected with the parameters A^F^) by the relations 

A' 2 (F^) = -^A| 3 , (4) 

where Cf = 4/3 and (3q is the first coefficient of the QCD /3-function. Both ways are 
completely equivalent if the number of active fermions in the expression for flo does not 
depend on Q. Meanwhile in order to provide self-consistency of the analysis, we changed n/ 
in Eq. from 3 to 4 at Q equal to the c-quark mass m c = 1.5 GeV. For this reason the 
value of A' 2 depends on Q in our case, although the numerical effect is inessential. The value 
of a s {Mz) was fixed at 0.118, which is close to the world average of Ref. ||. As one can see 
in Fig.l the analyzed data are insensitive to the parameter A 2 and we fixed it at the value 
of 1 GeV 2 inspired by the results of Ref. |I| on the analysis of charged leptons DIS data. 
The systematic errors on the data were accounted for in the covariance matrix approach, 



described in Ref. Vt 



Table 1: The results of the fit of the IRR model to the data from different neutrino experi- 
ments. The value of x 2 over the number of data points (NDP) is given in the last column. 



Experiment 


A 2 [GeV] 


A 2 , [GeV] 


X 2 /NDP 


IHEP-JINR 


0.69 ±0.37 


1. 


3/12 


CCFR 


0.36 ±0.22 


0.91 ±0.77 


253/222 



The results of the fit are given in Table 1 and in Fig. 1. The obtained contribution to 
X p(y +U ) N j s negative which supports the earlier observation of Refs. [fj, |3| and is in agreement 
with the results of Refs. || ||. The HT contribution to F^ U+U ^ N is negligible in the whole 
region of x spanned by the data. The value of A| is determined from our data with the 
50% accuracy. For the comparison, in the NLO QCD fit to the CCFR data on the structure 
function xFf +v)N the value A' 2 {F» N ) = -0.12 ± 0.05 GeV 2 was obtained in Ref. §. This 
estimate did not account for the systematic errors in the data, while the estimate accounting 
for systematics is A' 2 (F^ N ) = —0.10 ± 0.09 GeV 2 0, i.e., our result is the most precise 
estimate of the IRR normalization parameter at the moment. 

The change of A 2 under variation of a s (Mz) by ±0.003 is 0.055 GeV 2 and we consider 
this shift as a theoretical error in A 2 . Another source of the theoretical error comes from 
the uncertainty due to the effect of the higher-order (HO) QCD corrections to the LT term. 
These effects may be especially important for our study since the data at rather low Q are 
involved in the analysis. The HO corrections generally make the Q dependence of the LT 
contribution steeper, and correspondingly lead to the decrease of the HT contribution. In 
order to estimate the uncertainty due to neglected HO corrections, we repeated the fit with 
the QCD renormalization scale changed from the nominal value of Q to 2Q (see Ref. [pi| for a 
detailed argumentation of this approach). The obtained shift in the value of A 2 is 0.15 GeV 2 . 
Note that the significant part of the error in the world average of a s (Mz) also come from the 
uncertainty due to neglected HO QCD corrections. For this reason the two considered sources 
of theoretical errors are correlated. Having no possibility to account for this correlation, we 
just combine both errors in quadrature and estimate the total theoretical error in A3 as 
0.16 GeV 2 . One can see that the uncertainty in A 3 is dominated by the experimental error, 
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Figure 1: The x dependence of the measured structure functions F 2 (upper) and 

X p( u + U ) N (lower). The average values of Q 2 (GeV 2 ) for the x bins are given in the up- 
per plot. The full curves give the result of the LT+HT fit to the data, the dashed curves 
correspond to the la bands of the HT contributions obtained from the fit. 
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moreover accounting for the correlations of the separate sources of the theoretical error would 
decrease the latter. 

3. We also extracted the HT contribution to F 2 % from the CCFR data of Ref. [Q using 
the approach with a s fixed. The value of the parameter A3 obtained from the NLO QCD 
analysis of those data with x < 0.7 is given in Table 1. The CCFR data are sensitive to 
the parameter A 2 too, although the precision is poor and the fitted value is comparable with 
zero within the errors^. The theoretical error in A| estimated in the same way as for the 
analysis of our data is 0.12 GeV 2 . Results for both experiments are comparable within the 
errors and combining them we obtain the average 

A3 = 0.44 ± 0.19 (exp) GeV 2 . (5) 

In order to check universality of the IRR model scales with respect to the specific choice 



of structure function in the DIS process, we compared this value with the results of Ref. [16| 



on the analysis of the charged-leptons DIS data. Since in Ref. |16|] the results are given in 



terms of parameter A' 2 , we transformed our average fl5|) using Eq. (f|). As a result, we obtain 
that for rif = 3 

A' 2 {F£ N ) = -0.130 ± 0.056 (exp) GeV 2 . (6) 

This value is smaller, than A' 2 (F^ N ) = -0.2 GeV 2 , given in Ref. [|1(J, although within 
the errors both values are comparable. More precise conclusion about universality of the 
IRR model scales may be derived from the analysis of experimental data with improved 
statistics, which have been collected using the IHEP-JINR Neutrino Detector with a different 
configuration of the neutrino beam channel. The analysis of these data in order to extract 
the structure functions is currently in progress. The data from the NuTeV collaboration 
2(J, after their processing have been completed, may also be used to improve the precision 



of the IRR scales determination. In far sight a potential neutrino factory would allow for a 
detailed cross-check of the IRR model predictions and, in particular, the determination of 
the IRR scales with an accuracy of several percent 



In conclusion, we extract the high-twist contribution to the neutrino-nucleon structure 
function xF^ u+ ^ N from the analysis of the data collected in the first runs of the IHEP-JINR 
Neutrino Detector at the IHEP U-70 accelerator. We observe the negative HT contribution to 
the structure function xF^ u+u ^ N which supports the earlier observations. The normalization 
scale of the IRR model extracted from the combined analysis of the IHEP-JINR and CCFR 
experiments is about la lower than the one extracted from the data on the structure function 
F 2 N for the DIS of charged leptons. 
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